Theoretical study of the (3 x 2) reconstruction of /3-SiC(001) 
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By means of ab initio molecular dynamics and band struc- 
ture calculations, as well as using calculated STM images, 
we have singled out one structural model for the (3 x 2) re- 
construction of the Si-terminated (001) surface of cubic SiC, 
amongst several proposed in the literature. This is an alter- 
nate dimer-row model, with an excess Si coverage of | , yield- 
ing STM images in good accord with recent measurements 
[F.Semond et al. Phys. Rev. Lett. 77, 2013 (1996)]. 
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The reconstructions of SiC (001) surfaces have been 
widely studied in the last ten years!, the characteriza- 
tion and understanding of growth mechanisms on the 
(001) substrate being prerequisites for technological ap- 
plications. In the case of Si-terminated surfaces, sev- 
eral reconstructions have been found to occur; (2 x 1), 
c(4 x 2) and (n X. 2j)n=3,5,7... periodicities haw been ob- 
served in LEEOrM, RHEEDB and STMETtll measure- 
ments. Both (2x1) andc(4x2) reconstructions pertain to 
a complete Si monolayer at the top (9si = l^rfs clearly 
indicated by all available experimental dataoO. Unlike 
those of Si(001), the reconstructions of SiC(OOl) axe-char- 
acterized by weakly bonded, flat dimers ((2xl)E3E3) or 
by alternating symmetric dimers with different heights 
(c(4x2)E3c£l). Adsorption of additional Si produces suc- 
cessive (nx2) reconstructions as a function of Si coverage, 
including (7x2), (5 x 2j ,1-aJ^.a.combination of (5 x 2) 
and (3 x 2) periodicitiesffTfflEM 

The (3 x 2) reconstruction seems to be the last stage 
before self-limitation of growthEJ. Its atomic config- 
uration and electronic structure are not clearly estab- 
lished, though they have been intensively investigated. 
Three different atomic configurations, depicted on Fig. 1, 
have been suggested in the literature. In tbc Double 
Dimer-Row (DDR) model, proposed by DayanB-apd-Bas.- 
sibly supported by other experimental studieaJOOo, 
there are two Si ad-dimers on top of the full Si layer 
(Fig. l.a). The resulting coverage dsi — § is in con- 
tradiction with the measured 9si value of i reported 



by several groupsQ'Eit3. The straightforward extension 
of this model to the (5 x 2) reconstruction is also in- 
consistent with the measured coverages. Moreover, this 
model is not supported by some STM studieaffl. The 
DDR is favored by empirical molecular dynamics (MD) 
simulations^, but the Tersoff potential used in these cal- 
culations is known to give a poor description of /3-SiC 
surface reconstructions^. Another model, the ADDcd 
dimer-row (ADD), was first suggested in an early study 
by Hara et ajj. This configuration, with one Si ad-dimer 
per unit cell (Fig. Lb), corresponds to the measured cov- 
erage for the (3 x 2) and (5 x 2) reconstructions. However, 
though it appears ceasistent with several experimental 
data, both empiricalEil and ab initios calculations have 
shown that it is not energetically favored. Furthermore, 
STM investigations do not support this modelErQ. An- 
other i coverage model, the ALTernate dimer-row (ALT) 
(Fig. l.c), was proposed by Yan et aOE3j_.ILhis con- 
figuration is .supported both by calculationsErEa and by 
STM studiesQft However, it cannot account for the ob- 
served relation between single domain-. LEED patterns 
with (2 x 1) and (3 x 2) periodicitieaMl. It also fails 
to explain the (3.-X 1) reconstruction observed after O 
or H adsorptioncrO. Note that all three models involve 
Si ad-dimers which are perpendicular to the dimers on 
the underlying Si surface. Indeed, previous calculations 
have shown that a single parallel ad-dimer is energetically 
much less favored than a perpendicular oneE3. 

All three models show some discrepancies either with 
existing experiments or calculations, but none of them 
can be safely ruled out, owing to the lack of consistency 
between all available data. In this contribution, we re- 
port the results of self-consistent ab initio total energy 
calculations for all of the three structural models, includ- 
ing full geometrical optimizations. The surface energies 
are compared using Grand Canonical Potentials. The 
computed dispersion of electronic states, as well as STM 
images, are compared to experiments. Considering all 
the evidence, we conclude that the (3 x 2) reconstruction 
of the Si-terminated surface of SiC(OOl) is best described 
by the ALT model. 
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FIG. 1. Ball-and-stick representation of the relaxed atomic 
structures for the DDR (a), ADD (b), and ALT (c) mod- 
els, d and 5z are the distance and height difference between 
adatoms, respectively. Only the ad-layer and the first under- 
lying Si layer are shown for clarity. Bonds are drawn if the 
distance between atoms is smaller than 2.7 A apart. 

Our calculations were performed at T = within the 
Local Density Approximation, using ab initio molecu- 
lar dynamics codes employed in previous studies of SiC 
surfacesEj. Fully nonlocal norm-conserving pseudopoten- 
tials were used for Si (s and p nonlocality) and C (s 
nonlocality)Eil. The system was simulated by a slab in a 
periodically repeated supercell. The bottom layers were 
frozen in the p(2 x 1) configuration determined earlierEj. 
Two different sets of calculations have been performed. 
To determine the relaxed surface atomic structures, we 
used a (6 x 4) supercell with 8 atomic layers and a 10 A 
vacuum region (the total number of atoms is 176 (184) 
for ALT and ADD (DDR) models). The plane- wave en- 
ergy cutoffs for the wave functions and the charge density 
were 36 Ry and 130 Ry, respectively. Sums over occupied 
states were performed at the T point, which corresponds 
to 4 incquivalcnt k-points in the Brillouin Zone (BZ) for 
a (3 x 2) cell. Next, the electronic band structure was 
computed in a (3 x 2) supercell with 12 atomic layers 
and a 6 A vacuum region, using an extension of ab-initio 
MD codes to finite wave functions vectors, and keeping 
all atoms fixed. Atomic positions in the six top layers 
and the two bottom layers were taken from the preced- 
ing ab-initio MD calculations, whereas those in the four 
central layers were assumed bulk-like. Wave functions 



and charge densities were expanded in plane waves with 
cutoffs of 40 Ry and 160 Ry, respectively. In these calcu- 
lations, the electronic charge density was computed using 
8 special k-points in the.-BZ, generated according to the 
Monkhorst-Pack schemeEa. 

The relaxed atomic structures for the three models 
are shown on Fig. 1. In the DDR geometry, one ad- 
dimer is strongly tilted (Sz — 0.62 A) and has a short 
bond length (d = 2.26 A) while the other, weakly bound 
(d = 2.66 A), is almost flat (Sz = 0.03 A). The in- 
equivalence fifnthe two ad-dimers disagrees with simple 
expectationsa'til and with previous calculations by Kita- 
batake-yt al, who found two flat ad-dimers for the DDR 
modeO. Their use of the Tersoff potential could explain 
this disagreement, owing to the neglect of charge trans- 
fer between Si and C atoms. A single flat and weakly 
bonded ad-dimer (d — 2.62 A) is obtained in the ADD 
model, the geometry being close to that previously ob- 
tained by Yan et al in a calculation similar to oursE3. 
Note that two slightly different configurations are possi- 
ble within the same model, since the weakly bonded Si 
dimers in the underlying layer can be arranged e ither all 
on one side, as originally proposed by Yan et alcilO, or 
in a staggered pattern (see Fig. 1). Starting from dif- 
ferent configurations, out-calculations always converged 
to the staggered patterned. Finally, in the ALT model, 
the ad-dimer is strongly tilted (Sz = 0.5 A) and strongly 
bound (d =— 2.24 A), in good agreement with previous 
calculationsE3. The length of the weak Si dimers in the 
underlying surface layer is close to the value computed 
for the (2x1) reconstruction using the same method and 
a 4 x 4 supercellEj. 

In order to determine the most stable configuration, 
we have compared total energies. The ALT model is 
lower in energy than the ADD model by about 0.5 eV 
per (3 x 2) cell; this energy difference is clearly in fa- 
vor of the ALT, our error bar being 0.3 eV, as estimated 
from cutoff and BZ sampling tests. The drastic reduc- 
tion from the 3.6 eV energy difference quoted by Yan et 
alE3 is likely due to their poorer BZ sampling and, to a 
lesser extent, to the additional relaxation leading to the 
staggered pattern. A direct comparison with the DDR 
model is not possible since the corresponding Si coverage 
is different. This difficulty can be overcome by using the 
grand canonical schemeo. The computed surface energy 
differences as a function of the Si chemical potential fj,si 
are shown in Fig. 2. The value of /J,si for bulk silicon was 
calculated with an energy cutoff of 40 Ry and 32 special 
k-points in the BZ, whereas we used the heat of foxma- 
tion of SiC AH = 0.75 eV from a recent calculations, to 
obtain the Si chemical potential under C-rich conditions. 
Since (3 x 2) growth is likely to occur under Si-rich con- 
ditions, a precise determination of ^si is required only in 
that limit. Our results show that the ALT model is the 
most stable configuration over the entire allowed range of 
Si chemical potential. However, the energy difference be- 
tween ALT and DDR obtained under Si-rich conditions 
is only 77 meV, i.e. within our error bar. Consequently, 



the DDR model can not be definitely ruled out solely on 
the basis of total energy comparisons. 
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FIG. 2. Total energy difference per (3 x 2) unit cell as a 
function of the Si chemical potential fisi- 
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FIG. 3. Calculated constant-current STM images (bias 

V = -1 V) for the DDR (a), ADD (b), and ALT (c) models. 

Correlation with the atomic structure is explicitly shown via 

the superposed ball-and-stick representations of the atoms in 

the first surface layer (see fig. 1). 

Several experimental STM st«udjes.pf the (3 x 2) recon- 
struction are currently availableBl3'E3'E3. In order to com- 
pare the three different models, we have calculated filled 
states constant-current STM images within the Tersoff- 
Hamann approximationEU. Representative images are 
shown in Fig. 3. In both the DDR and ADD mod- 



els we find strings of peanut-shaped spots, originating 
from a slight overlap between maxima on adjacent flat 
ad-dimers. For the DDR model, additional maxima are 
located on the up adatoms of the tilted ad-dimers. The 
resulting images are incompatible with the experimental 
observations of a single oval spot stretched in the [110] 
direction per 3x2 cell. On the other hand, in the ALT 
model the spread out stretched spots located above up 
adatoms of the tilted ad-dimers are in accord with exper- 
imental STM images of filled stateso. 

Additional insight can be obtained from analysis of 
the electronic states within a few eV of the Fermi level. 
All photoemission measurements agree about the pres- 
ence of two occupied surface— states in the band gap, 
1 eV apart from each otheniSEiffiij. However, uncertain- 
ties exist about the location of these states with respect 
to the Valence Band Maximum (VBM). Recent angle- 
resolved photoemission spectroscopy (ARPES) measure- 
ments have shown that the dispersion of all identified, 
surface states, is very small (< 0.2 eV) along the [Tl0]lia 
and [110JI3O directions. Only the surface states of the 
DDR and ALT models have been considered here, the 
ADD model being higher in energy than the ALT model 
and exhibiting STM images which do not agree with ex- 
periment. 

We find that in the DDR model the surface is metal- 
lic, within the Local Density Approximation. The high- 
est occupied state, about 1 eV above the VBM at T, 
is mainly localized on the flat ad-dimer and has a 7r* 
character with respect to the dimer axis. Its dispersion 
is very small along the [110] direction (< 0.1 eV), but 
rather strong along the [TlO] direction (~ 1 eV). In the 
DDR model, we also find three additional surface states 
with energies between the highest occupied state and the 
VBM. Only one of them is localized above the up adatom 
of the tilted ad-dimer, and is essentially dispersionless. 
The other two states originate from backbond and dimer 
states of the underlying surface, and show strong dis- 
persions. The presence of dispersive states in the band 
gap is in disagreement with ARPES evidence and points 
against the DDR model. 

In the ALT model, the surface is semiconducting, with 
a direct gap at T of about 0.5 eV. The highest occupied 
state, 0.8 eV above the VBM at T, is localized on the up 
adatom of the tilted ad-dimer and has a strong 's' char- 
acter. This surface state has a small dispersion along 
both [110] and [TlO] directions (< 0.1 eV). Close to the 
VBM we find another state, lying 0.7 eV below the high- 
est occupied orbital. It is a it* state localized on the Si-Si 
dimer of the underlying surface which are not bonded to 
ad-dimers, and is nearly dispersionless (< 0.2 eV). This 
state is only present in the [110] direction. Except for the 
energy difference between the two highest surface states 
(0.7 eV vs. 1 eV), agreement with ARPES experiments 
is definitely better for the ALT than the DDR model. 

A discrepancy exists about the numbec. r ajid location 
ofj-occupied surface states or resonancest§Ei Yeom et 
allia argued that there should be a total of four states 



with pronounced surface character in their proposed 
DDR geometry. However, they considered the ideal non- 
relaxed geometry, with flat ad-dimers only, and, more 
importantly, they ignored backbond states and dimcr-likc 
states on the underlying surface which are also expected 
to exhibit surface character. These additional resonant 
states could be difficult to resolve because some are close 
in energy, and might have weak photoemission intensi- 
ties. In our calculation we could identify resonant states, 
however we did not attempt to systematically analyse 
their character and dispersion. 

Turning to unoccupied states, we find only one sur- 
face state in the band gap for the DDR model. It lies 
2 eV above the VBM at T, is localized around the down 
adatom of the tilted ad-dimer and has a predominant 
'Px' character. Its dispersion is about 0.1 cV (0.4 eV) 
along the [1 10] ([110]) direction. For the ALT model, 
two empty surface states have been identified. The low- 
est one, with energy 1.2 cV above the VBM at T, is a a- 
like state on the lone Si dimer of the underlying surface. 
It disperses about 0.1 eV along [110] and 0.7 eV along 
[110]. The other one, 1.8 eV above the VBM, is a 'p z '-like 
state on the down adatom of the tilted ad-dimer and has 
a very weak dispersion along both directions (< 0.2 eV). 

The combination of all our results indicates that the 
ALT model is the most suitable candidate, since it ex- 
plains the large majority of available measurements. The 
ADD and DDR models produce incorrect STM images. 
The ADD is energetically unfavorable, while the disper- 
sion of the surface states calculated for the DDR model 
is not compatible with ARPES measurements. 

To summarize, we have performed plane-wave pseu- 
dopotcntial calculations for three structural models of the 
(3 x 2) reconstructed /3-SiC(001) surface. In particular, 
relaxed atomic structures, surface energies, STM images, 
surface states and their dispersion have been calculated 
and compared with experiments. Our results strongly 
favor the ALT model and exclude the DDR and ADD 
models, although some ambiguities remain. More defini- 
tive conclusions could come from additional experimental 
studies, in particular investigations of unoccupied elec- 
tronic states, and a convincing confirmation of the Si 
coverage corresponding to the (3 x 2) reconstruction. 
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